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We transpose an idea of 't Hooft from its context of Yang and Mills' theory of strongly interacting 
quarks to that of strongly correlated electrons in transition metal oxides and show that a Hubbard 
model of N interacting electron species reduces, to leading order in N, to a sum of almost planar 
diagrams. The resulting generating functional and integral equations are very similar to those of 
the FLEX approximation of Bickers and Scalapino. This adds the Hubbard model at large N to the 
list of solvable models of strongly correlated electrons. 
PACS Numbers: 71.27.+a 71.10.-w 71.10.Fd 



is by itself reason enough for studying 



I. INTRODUCTION 

Recent numerical [|J and analytical j^] studies have shown that an extrapolation of the Hubbard model from two to 
N distinct electron species, and with its SU(2) invariance replaced by SU(N), conserves the nontrivial character of 
this model and this provides a motivation for further study of this extrapolation. Actually, the mere fact that SU(N) 
antiferromagnets have spin dimer or "RVB" ground states at N = oo 
the SU (N) extrapolation of the Hubbard model. 

Existing methods of treating such large N limits are of two kinds. The first is a quasi classical saddle point method 
used in the context of nonlinear sigma models || and slave operator constraints || and that was applied already to 
the SU(N) Hubbard model |Q, but without conclusive results, in the opinion of this author. The second approach is 
't Hooft's topological expansion || as applied to a SU (N) deformation of the Yang-Mills theory of strong interactions. 
In this latter approach, N fermions interact with N 2 bosons and a quasi classical or saddle point interpretation is 
therefore impossible. 

't Hooft's topological large A" expansion is less well known, but more profound than the large N saddle point 
method. In this paper we will use it to show that the SU(N) Hubbard model is solvable, for A~ 3> 1, in terms of 
planar diagrams that can be summed via coupled integral equations. 

Remarkably, our generating functional and associated integral equations turn out be very similar to those of the 
" fluctuation exchange approximation" (FLEX) of Bickers and Scalapino || . 

II. A TOPOLOGICAL CLASSIFICATION OF DIAGRAMS 

We consider the Hubbard model [jl0| that epitomizes the problem of strongly correlated electrons (see [jll) for a 
review) and extend it from 2 to A" species of electrons H : 



x,y x \a=l..N / 



(1) 



We have replaced U — > jf to have a stable limit as N varies and to regain the conventional meaning of U at A" = 2. 
A suitable chemical potential is implicit in all our expressions. 

A satisfactory treatment of the interaction between electrons and their own spin and density fluctuations is at the 
heart of the correlation problem. In the model of Hubbard (and in essentially any electronic system with finite range 
two body potentials) this interaction can be made explicit by means of a Hubbard-Stratonovich transformation: 

fj ip f3 + h.c. I dr (2) 



J D<j> a f) CXp - J + y^C^a 

const(T)exp-^- I ( V" ty*a$aV$>p + density ] dr 



Here "density" stands for terms ~ ip^a that can be absorbed in a redefinition of the chemical potential and which 
we may therefore ignore. To avoid a proliferation of unphysical fields, we impose 0* g = <f)p a . We then find the 
following bare correlators for <p a g and ip a 



1 



< 0a/3(l)0Q'/3'(2) >o= ^6 aa '5pp'6(l, 2) and < Q| 30 a '/3' >0= ^a/3'<5/3a'<Kl, 2) (3) 

< Va(i)^( 2 ) >°= <w> (stt^;) (1,2) 

where 5(1, 2) is a shorthand for 8{t\ ~ T i)S-^f i -jf 2 - With the help of the spin and density fluctuations 4> a p the partition 
function of the repulsive Hubbard model can be written as 

Z = J Dq^Dyj exp - / ' dt [ £ (C^V-a + ^a/j) + H j (4) 

V a,f3,x J 

From eqs(|^^) the perturbative diagrams associated with the above Hamiltonian may be drawn in terms of single line 
propagators for electrons, double line propagators for the bosons and a three point vertex 2^/^"^* (ftapipfl at which 
a particle hole excitation splits into its constituents, see figure(l). 

/3 
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FIG. 1. Bare propagators and bare vertex in the SU(N) Hubbard model 

Using an argument that 't Hooft originally applied to an SU(N) deformation of Quantum Chromodynamics [§J, 
we now show that any perturbative diagram made up of such propagators and vertices carries a weight N x where 
X = 5*2 — Si + Sq is the Euler characteristic |Q of its associated topological surface. Here Sq stands for the number of 
vertices, Si for the number of propagators and S% for the number of "area like" pieces associated with each diagram. 
To define an "area like" piece, we follow an index line until it closes. The topological surface associated with each 
diagram then consists of a collection of loops that are glued together along segments that belong to double lines 
representing boson propagators. In fact, this is precisely the type of representation of a surface that is used in 
elementary topology |l2|] where for instance a torus is represented in terms of a rectangle with identified opposite 
sides. 

By definition then, each closed index loop of a given diagram contributes one unit to 5*2. We may count the 
endpoints of each propagator either from the propagator point of view (25*1) or from the vertex point of view (3<So) to 

obtain the topological relation Si — §<So- With a factor ^J~^ from each vertex and with a factor N from each closed 
index loop, a diagram carries a weight factor of 

weight = U^N S ^ 5l = |5 ° U%N* (5) 

The Euler characteristic of an orientable closed surfaces with h handles and o openings is x — 2 — 2h — o. The 
last equation implies that the "planar" diagrams that can be drawn in a plane without intersections of lines are the 
important ones, while extra openings and handles are suppressed by factors of -i and to, respectively. 

Although it is difficult to intuitively understand the intrinsic topology of the topological surface associated with 
any given diagram, it is easy to draw a surface with a minimal number of handles and holes on which a given diagram 
can be drawn. Its intrinsic geometry and Euler characteristic coincide with that of the topological surface defined via 
single and double line propagators and vertices. We illustrate this with a nonleading contribution to the electron self 
energy in figure (2) where the red lines represent boson propagators and the black line an electron. 
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FIG. 2. Nonplanar and therefore nonleading contribution to the fermion self energy 

The extra handle is needed here for the boson propagators not to cross each other. From the factor N x of eq(|5|) 
and the expression for the Euler characteristic we see that this diagram is suppressed, relative to the planar ones 
(and without crossings of boson propagators) by a factor of to. Of course we arrive at the same conclusion also by 
carefully drawing this diagram in terms of lines and double lines and by counting powers of N from vertices and index 
loops. 

Finally, we must also mention a difficulty in the present SU (N) extrapolation of the Hubbard model. The diagram 
of figure (3) without the large half circle contributes to the bosonic self energy and gives rise to scalar or trace 
contribution ~ a>3 N i' ' in the < <i> a fi4 l * a i 31 > propagator. 




FIG. 3. Self energy correction due to a scalar part in the boson propagator 

This means that the tensor character of the full boson propagator is different from that of the bare one in eq(|^). 
Luckily, this term is suppressed at large N, as we shall see. The small circles in figure (3) indicate vertices, of which 

there are six (for the electron self energy) with a factor of ^J~^ each and hence this diagram gives a contribution of 

order O(jks) to the electron self energy. We shall, therefore, ignore the scalar or trace part of the < (fxf)* > propagator 
on the grounds that (i) its effects are small at large N and (ii) because it represents only uninteresting charge density 
fluctuations JDi]. 



III. COUPLED INTEGRAL EQUATIONS AND GENERATING FUNCTIONAL 

Expanding the cubic interaction in eq(^J) to second order (or by appealing to standard diagram lore explained in 
textbooks jnj) we find the following lowest order electronic self energy 

4U 4J7 

Eaxa, (1,2) = "jy < ¥Va 3 (l)V>a 3 (l) '^(2)^(2) >= — £> QlQ3 | QlQa (l, 2)G Q3Ql (l, 2) (6) 

G Q3Q4 (1,2) = < T {ip a3 (l)ip* 4 (2)} > and -D QlCt3 | Q4Q2 (1, 2) =< T {ip 

Q1Q3 V x 7 Vet 4 a 2 (2)}> 

We promote this expression to the complete skeleton one by using full propagators in it. This, however, ignores vertex 
corrections. The first vertex correction to the electronic self energy E is suppressed by a factor of i and may be 
ignored at leading order in N. 

A more graphic way of arriving at eq(|^) is by considering the sum of planar or " rainbow like" contributions to the 
electronic self energy and by removing the outermost scalar propagator or rainbow. 
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To exploit eq(^) we use the tensor structure of D and G which is (ignoring the scalar or trace contribution to D) 

-^(1) 2) aiQ3 | Q4a2 = (5 aiQ2 ^ CK3a4 -D(l, 2) and 11(1, 2) aia3 | Q4Q2 = 6 aia2 5 a3a4 Tl(l 7 2) (7) 
Ga 3 a 4 (l, 2) = S aaa4 G(l, 2) and £ aiQ2 (1, 2) = 2(1, 2)<5q, iQ2 

and we find 

£(1,2) = 4[/ J D(l,2)G(l,2) (8) 

Starting again from second order perturbation theory for II and promoting it to a skeleton diagram (which is correct 
to leading order in N) we find a boson self energy of 

4U 4U 
n ai/W2 (l,2) = — < ^(1)^(1)^(2)^(2) >=-— <J Ql/3 A 2/3l G(l,2)G(2,l) (9) 

n(l,2) = -^G(l,2)G(2,l) 

where the minus sign reflects the fermion loop in this diagram. To be sure of the self consistency of the integral 
equations for £ and II and, at a later stage, for thermodynamics and linear response, we need a functional that 
generates the equations for £ and II. To find it, we consider the lowest order diagram for the logarithm of the 
partition function and promote it also to a skeleton diagram: 

f = ¥ E<< (!)*««<« (!)V>a 2 (1)C 3 (2)^4 (2)^ 4 (2) > (10) 
2[j 

= ^ t D aiot2 \u3a4 (1; 2)G a2C[3 (1, 2)G aiai (2,1) 

Here ^2 stands for a sum over repeated indices and for integration over Matsubara time and the extra factor 1/2 is 
due to the symmetry of this diagram. Combining eqs(|7|,|i~0|) one recognizes that F contains two closed index loops and 
one obtains 

F=-2C/^D(1,2)G(1,2)G(2,1) (11) 

with full Green's functions D, G. By general field theory or many body lore []l5f our previously obtained equations 
for £, II should come out by varying F with respect to D, G: 

E(l) 2) = -tJtu E -WD(1', 2')G(1', 2')G(2', l') (12) 
n(1 ' 2) = ^wkl) E -WD(1', 2')G(1', 2')G(2', l') 

(diagrammatically speaking, functional differentiation applied to vacuum diagrams opens propagator lines in all 
possible ways). Dyson's equations G^ 1 = G^ 1 — £ and D^ 1 = D^ 1 — II for the self energy may also be written in 
differential form: 

S(l, 2) = -^A^Tr (log G-G/G ) (13) 

Taken together, eqs([l2],[l3]) permit the introduction of a new functional f2 that is stationary with respect to variations 
in G and D: 

n = Tr (log G - G/G ) - ^Tr (log D - D/D ) + £ 2UD(1, 2)G(1, 2)G(2, 1) (14) 

1,2 

sn _ sn 

5G~5D~ 

That Q represents the logarithm of the partition function of the system under consideration can be proved by general- 
izing the arguments of Luttinger and Ward |lf| to coupled bosons and fermions. The factor N/2 reflects the presence, 
in our system, of ./V 2 real bosons as compared to only N complex fermions. 



4 



IV. THE PHYSICAL MEANING OF THE INTEGRAL EQUATIONS 



To gain a first understanding of the integral equations one may solve them iteratively in powers of 1/N. From 
eqs(|| ||) we see that one may set G = Go to leading order and then one finds 

G = G + O(^) (15) 

11(1, 2) = -^G (l, 2)G (2, 1) + 0(^_) 

ft is now clear that II represents paramagnon type spin fluctuations (see |l6| for an appraisal of such theories and 
for further references) and we must compare with the FLEX approximation of Bickers and Scalapino |^] that was 
designed to deal with particle hole fluctuations in correlated systems and, in particular, with spin fluctuations in the 
Hubbard model. 

Since the FLEX approximation is stated most succinctly in terms of a free energy or generating functional we 
compare with flT] ] where the FLEX approximation to the free energy of the Hubbard model is written down explicitly. 
The first thing we notice from jlTj is that the FLEX functional involves only G but not D. To compare our functional 
with the FLEX one, we must eliminate D from it by using the saddle point condition eq(||). Together with Dyson's 
equation for D, this leads to a drastic simplification of the functional Q : 

N 

[n] sa =Q = const + Tr (log G - G/G ) - — Tr (log D) (16) 

where D^ 1 = Dq 1 — U. Readers who are suspicious of the simplicity of this functional may be reassured by the fact 
that a functional that includes the first subleading order (and which is not given here) has a more subtle mathematical 
structure. 

To compare the functional of the last equation with that of FLEX we note that according to |l7j the FLEX 
functional contains no linear term in n, while ours does (when expanding in n or 1/N) so the two functionals do not 
coincide. Another difference is that FLEX violates crossing symmetry |Q while the large N approach does not. 

Such differences between the two approaches are not surprising in view of the fact that the Kadanoff-Baym pre- 
scription for writing down conserving approximations is non unique ]T^ ] while there is no ambiguity in the large N 
approach as described here, except that one may mix spin and charge fluctuations with an arbitrary mixing angle. 
This mixing angle ceases to be arbitrary only if any one of the fluctuating fields develops an expectation value. How- 
ever, because we do not expect any condensation of local pairs ip a (x)ipp(x) in the repulsive Hubbard model, there is 
no reason to introduce such a mixing of channels in the present context. 



V. CONCLUSIONS 



Our main result is the observation that the SU (N) extension of the Hubbard model can be solved, at large TV, 
by a transcription of 't Hooft's planar diagram technique from SU(N) Yang-Mills theory to this model. At leading 
order in N the functional f2 of cq(|l4|) must be minimized via the saddle point equations (8,9) for E and n. Because 
our generating functional and integral equations are quite similar to the FLEX equations of Bickers and Scalapino we 
expect that their solutions should be qualitatively the same. 

Since the ground state of the undoped SU (N) theory for large TV is a collection of RVB type dimers our arguments 
constitute a bridge between Anderson's early RVB ideas [Q and the FLEX method. 

One of the shortcomings of the SU(N) extrapolation of the Hubbard model presented here is that it is unlikely to 
describe antiferromagnetic order, such order being unnatural for more than two electron species. It is probably no 
coincidence that the absence of antiferromagnetic order and of a spin gap is also the main shortcoming of the FLEX 
approximation. 

Clearly, the approach presented here must be improved in a way that permits a smoother extrapolation back to 
the electrons of our physical world with their two spin orientations and their (anti) ferromagnetic correlations. 
Acknowledgements 

I am indebted to P. Schuck and J.-M. Robin for useful correspondence and to A. Hewson, J. Keller, T. Pruschke 
and M. Zoelfl for informative discussions. I thank T. Dahm for extensive correspondence on the FLEX approach and 
I am indebted to S. Villain- Guillot for numerous discussions of the manuscript. Finally I wish to acknowledge the 
marvelous inspiration provided by the 1999 Hvar School on Correlated Electrons. 



5 



[1] O. Gunnarsson, Z. Physik B104 (1997) 279. 

[2] R. Assaraf, P. Azaria, M. Caffarel, P. Lecheminant, Phys. Rev. B 60 (1999) 2299. 

[3] N. Read and S. Sachdev, Nucl. Phys. B315 (1989) 609; D. S. Rokshar, Phys. Rev. B42 (1990) 2526 and (for O(oo) instead 

of S77(oo)) D. Foerster and F. Triozon, Phys. Rev. B56 (1997) 8069. 
[4] For details and references on the RVB idea see P.W. Anderson, "The Theory of Superconductivity in the high T c Cuprates" 

Princeton University Press, Princeton, N.J. 1997. 
[5] H. E. Stanley, Phys. Rev. 176 (1968) 718. For examples of both versions of the large N approach, see the reprint volume 

"The large N Expansion in Quantum Statistical Physics", edited by E. Brezin and S. Wadia, World Scientific, Singapore 

1993. 

[6] for reviews see: A. C. Hewson, "The Kondo Problem to Heavy Fermions", Cambridge University Press (1993) and N. E. 

Bickers, Rev. Mod. Phys. 57 (1987) 845. 
[7] I. Affleck, B.J. Marston, Phys. Rev. B. 39 (1989); I am indebted to P. Lecheminant for calling this paper to my attention. 
[8] G. 't Hooft, Nucl. Phys. B72 (1974) 461 and Nucl. Phys. B75 (1974) 461. 
[9] N. E. Bickers and D. J. Scalapino, Ann. Phys. (N. Y.) 193 (1989) 206. 
[10] J. Hubbard, Proc. Roy. Soc London A 276 (1963) and A 281 (1964) 401. 

[11] P. Fulde, "Electron Correlations in Molecules and Solids", Springer Series in Solid-State Sciences 100, N.Y. 1991. 
[12] See for example J. Stillwell, "Classical Topology and Combinatorial Group Theory", Graduate Texts in Mathematics 72, 
Springer, N.Y. 1993. 

[13] There is a related complication with the vertex that changes its tensor character at higher order in 1/N but the corre- 
sponding contribution to the electronic self energy is suppressed by a factor of 1/N 2 

[14] A. A. Abrikosdov, L. P. Gorkov and I. E. Dzyaloshinski, "Methods of Quantum Field Theory in Statistical Mechanics, 
Dover, N. Y. 1997. 

[15] G. Baym, Phys. Rev. B 127 (1962) 1391; J. M. Luttinger and J. C. Ward, Phys. Rev. 118 (1960) 1417. 
[16] Y. M. Vilk and A.M. Tremblay, J. Phys. (France) I 7 (1997) 1309. 
[17] J. W. Serene and D. W. Hess, Phys. Rev. B 44 (1991) 3391. 
[18] N. E. Bickers and S. R. White, Phys. Rev. B 43 (1991) 8044. 



G 



